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Abstract. If X and Y are independent, Y and Z are independent, and so are X and Z, 
one might be tempted to conclude that X, Y, and Z are independent. But it has long 
been known in classical probability theory that, intuitive as it may seem, this is not 
true in general. In quantum mechanics one can ask whether analogous statistics can 
emerge for configurations of particles in cerfain fypes of entangled states. The explicit 
construction of such states, along with the specification of suifable sets of observables 
that have the purported statistical properties, is not entirely straightforward. We 
show that an example of such a configurafion arises in the case of an N-parficIe GHZ 
sfate, and we are able to identify a family of observables with the property that the 
associated measurement outcomes are independent for any choice of 2,3,..., N - 1 
of the particles, even though the measurement outcomes for all N particles are not 
independent. Although such states are highly entangled, the entanglement turns out 
to be 'fragile', i.e. the associated density matrix has the property that if one traces 
out the freedom associafed with even a single particle, the resulting reduced density 
matrix is separable. 
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1 . Introduction 

The notion of statistical dependence can at times be counterintuitive, leading to 
surprising conclusions even in classical probability lH]]. In the case of a quantum 
system consisting of two or more particles (or subsystems), the intrinsic quantum 
statistical dependence of the various constituents upon one another, otherwise known 
as 'entanglement', can admit features—such as violation of the Bell inequality— 
that do not have obvious classical analogues. Indeed, ever since Schrodinger's 
introduction of the idea of 'entanglement' ©[Si, and the concurrent work of Einstein, 
Podolsky and Rosen on its paradoxical aspects IIU, a great deal of effort has been 
made by researchers to understand the implications of this fundamental notion, and 
to seek applications (see, e.g., articles in [3 and references cited therein). 
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The purpose of the present paper is to point out the peculiar statistics arising from 
measurements of a certain set of observables when a system of particles is prepared in 
a particular type of entangled state. Our considerations lead to a quantum counterpart 
of a surprising result in classical probability due to Bohlmaim and Bernstein [Zl: to 
wit, that it is possible to construct a situation in which while X and Y are independent, 
Y and Z are independent, and X and Z are independent, nevertheless X, Y, and Z are 
not independent. We show, more generally, that when measurements are made 
on a quantum system consisting of N sptn-| particles, a version of the Bernstein 
phenomenon arises when the state of the system belongs to a particular class of fully- 
entangled states that we shall refer to as 'Bernstein states'. We construct this family 
of states explicitly and show that no other states give rise to Bernstein statistics. The 
elements of a special subclass of these states turn out to be equivalent, modulo local 
unitary transformations of the pure phase type, to N-particle GHZ states BBl. While 
such GHZ-equivalent states are known to be strongly entangled, they are at the same 
time 'fragile' in the sense that if one traces over the contribution of any one of the 
constituent particles then the entanglement is destroyed—that is to say, the resulting 
reduced density matrix of the remaining N - 1 particles is separable. 

The paper is organized as follows. In Section|2]we review the classical Bernstein 
distribution, and construct an example of a three-particle quantum Bernstein state by 
analogy with the classical case. This state is shown to be fragile. In Section |3l the 
three-particle example introduced in Section|2]is generalized by use of a combinatorial 
argument to the construction of a special N-particle Bernstein state for each N > 3. 
In Section HI we show by use of a basis transformation that the special N-particle 
Bernstein states are equivalent to a class of N-particle GHZ states and that these 
states are fragile. In Section|5l the N-particle Bernstein states constructed in Section|3] 
are shown to be unique in having the Bernstein distribution, up to the inclusion of a 
set of 2^“^ phase factors. In Section [6] we show that for N > 3 the space of Bernstein 
states that are equivalent modulo local phase transformations to the special Bernstein 
states, for which the phase factors take the value unity, is a proper subspace of the 
space of general Bernstein states. In the three-particle case, the special Bernstein states 
and the general Bernstein states are equivalent modulo local phase transformations. 
In Section [zl we construct a state with an inhomogeneous Bernstein distribution and 
show that the fragility property holds in the homogeneous limit. Finally, in Section |8| 
we show how the so-called Mermin paradox jUllIll can be extended by the use of our 
methods to the case of an N-particle GHZ state. 

2. Bernstein statistics 

For the benefit of readers less well acquainted with the relevant ideas, first we 
introduce the Bernstein distribution jjZlllQl. Gonsider a set of four cards, on which 
the following sequences of numbers are printed: (110), (101), (Oil), and (000). The 
four cards are placed in a bag, from which one card is selected at random. Let X 
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be the event that the first digit on the selected card is 1, let Y be the event that the 
second digit on the selected card is 1, and let Z be the event that the third digit on the 
selected card is 1. Clearly, the probability of X is the probability of Y is and the 
probability of Z is Similarly we find that the probability that X and Y occur is the 
probability that X and Z occur is j, and the probability that Y and Z occur is Thus, 
the triplet (X, Y, Z) is pairwise independent. Now, if X, Y and Z were independent, 
the probability of all three events occurring would be Since the probability of all 
three events occurring is zero, it follows that X, Y and Z are not independent. 

A quantum version of the Bernstein phenomenon can be constructed as follows. 
We consider a system of three spin-^ particles. Write Onx, Ony, (J„z for the spin 
operators associated with particle n, where n = 1,2,3, and write ti*^, 
for the corresponding projection operators onto the associated eigenstates. Thus, 
'^tx - 5(1 ± ^nx), Ti^y = 5(1 ± Onz), and 71*2 = ± ^nz)- In what follows, we work 

mainly with Onz eigenstates. Thus, for example, we write | tTT) for the state vector of 
a system consisting of three spin-l particles in Oiz, Oiz, Osz spin-up states. An explicit 
example of a Bernstein state in the case of a system of three spin-1 particles is 

|B™> = 5(ITTi> + ITa> + liTT> + l ). (1) 

Evidently, is neither a product state nor a partial product state, so it represents a 
totally entangled state. The structure of the expression on the right side of dH) allows 
one to verify the Bernstein property for spin measurements in the z direction. To this 
end, let us write (- 1 - • •) for the event that the outcome of a measurement of 0\z on the 
first particle is spin up; similarly, let us write (-1- -I- •) for the event that the outcomes 
of measurements of a\z and cr 2 z on the first and second particles, respectively, are both 
spin up, and so on for the various other possible outcomes for single, double, and 
triple spin measurements. One sees that when the system is in the state the 
probabilities of the outcomes of the various single spin measurements are given by 

P(+ . .) = P(. + .) = P(. • +) = 1 (2) 

along with similar results for spin down. Then for the probabilities of the outcomes 
of spin measurements on any pair of the three particles, we have 

P(++.)=P(+.+)=P(.++) = i, (3) 

with similar results if any 'up' is replaced with a 'down'. In particular, we find that 

P(+ -!-•) = P(-i- • •) X P(« -I- •), (4) 

and that similar decompositions into products of single spin measurement outcome 
probabilities hold for all possible outcomes of the three pairs. It follows that the 
outcomes for Uiz, and 032 are statistically independent for any pair of these 
measurements. On the other hand, if one measures the spins of all three particles the 
situation changes. In particular, we find that the probability of three 'ups' is given by 

1 

P(+ + +) = 0 ^ P(-i- • •) X P(« -I- •) X P(« •-!-) = -, 

8 


( 5 ) 
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which shows that the outcomes of the three spin measurements are not independent 
even though they are pairwise independent. In conclusion, one sees that when 
a quantum system is prepared in the state the statistics of the outcomes of 

measurements of the three spin operators Oiz, and Usz are those of Bernstein. 

The Bernstein state ([D is remarkable in another respect as well: the associated 
pure density matrix lies in the subclass of three-particle states having the 

property that if one traces out the degrees of freedom belonging to any one of the 
three particles then the resulting reduced density matrix is separable. This is the 
so-called 'fragility' property that we alluded to earlier. To verify that is 

fragile, we take the trace over, say, the first particle Pi, to obtain the reduced density 
matrix p 23 = trpi associated with the system consisting of the remaining 

two particles. Note that if, in the original three-particle state, one were to measure 
the first spin in the z direction, then conditional on an 'up' outcome the state of 
the remaining two particles would be -^ ( Iti.) -i- UT) )• Likewise, if the outcome of 
the measurement of the first spin is 'down', the state of the remaining two particles 
would be (ItT) + lii))- In both cases, the resulting two-particle state is maximally 
entangled according to the standard geometric measure of entanglement HTTl . On 
the other hand, if we trace out the first particle, the reduced density matrix of the 
two-particle system that results takes the form 



10 0 1 ' 
0 110 
0 110 
10 0 1 , 


( 6 ) 


when it is expressed in the basis ( ITT)/1 Ti)/UT)/1 i-i) )• The fact that this state is 
separable—in other words, that it can be represented as a mixture of disentangled 
states—can be verified by use of the Peres-Horodecki criterion 11121 [13]], which says 
that p 23 is separable if the partial transpose of p 23 has nonnegative eigenvalues. But 
it is clear that p 23 is invariant under partial transposition of either subsystem, and it 
has the noimegative eigenvalues j and 0, each with multiplicity two. It follows that 
the reduced density matrix is separable. Indeed, it is a remarkable fact that a mixture 
of two maximally entangled pure states can nevertheless be separable. This can be 
seen somewhat more directly by use of the identity 


P23 = 5(in> + iu>)(mi + aii) + 5(in> + iiT>)(mi + aTi) 

where |<—) = (jt) -i- ]],)) and j —>) = (|t) - jj.)), which expresses the reduced state 

p 23 as a mixture of entangled states and also as a mixture of disentangled states. 

Proceeding analogously, it should be evident then that similar results to those 
described above are obtained if we trace out the second particle P 2 or the third particle 
P 3 . From this, we conclude that the Bernstein state |B^^^) defined by ([1]) is indeed a 
fragile state in the sense that we have indicated. 
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A useful consequence of Q is that the independence of the measurement 
outcomes of cr2z and 02,2 can be seen to follow directly as a consequence of the 
symmetry of ([Z| with respect to the left and right arrow states. One can verify, 
for example, that tr P23 = tr p23 = 5, and tr p23 = which are the 

correct probability relations for independence of the indicator functions for the second 
particle having spin up and the third particle having spin up. Here we have used 
the fact that when measurements are made only on particles 2 and 3 , it suffices to 
make calculations using the reduced state P23. On the other hand, it should be noted 
that the separability of a two-particle density matrix such as p23 does not imply that 
the measurement outcomes of all bipartite observables are statistically independent. 
That can only happen when both the density matrix and the relevant observables 
factorize. In particular, the reduced state p23, which is separable, but not factorizable, 
gives rise to correlations. For example, we have tr 71+^ p23 = tr p23 = and 
tr 7X2^ 7X3^ p23 = \ ^ tr 7X2^ p23 tr n^^p23 = Hence, the outcomes of 02% and 03^ are 
positively correlated, but not the outcomes of a2z and 03^, which are independent. 

3. N-particle Bernstein states 

The example considered in the previous section can be generalized in various ways. 
First, we can ask whether the state dU is in some appropriate sense unique in admitting 
the Bernstein distribution for the indicated measurement outcomes. One possible 
generalization would be to consider a situation where the spin measurements on the 
various particles are taken with respect to a different choice of spin axis for each 
particle. But this set-up is equivalent to the original one in the sense that the outcome 
statistics are preserved if we make local unitary transformations on each of the spin 
operators in such a way as to line up the axes in the direction of a single choice of 
z-axis, and if at the same time we make the corresponding conjugate local unitary 
transformations on the various particles composing the overall state of the system. 
Thus without loss of generality we may assume, as we have done in the previous 
section, that the axes of the various spin operators are aligned along a fixed given 
axis, which we have taken to be the z-axis. 

On the other hand, even if we stick with the given set of spin measurements 
we note that the Bernstein distribution still arises if we generalize the state by 
inserting unimodular phase factors in front of each term appearing in ([T]). Before 
exploring the ramifications of the resulting degrees of freedom in this apparently 
larger family of three-particle Bernstein states, it will be convenient to consider a 
further generalization of the Bernstein state dD, namely, that arising when we consider 
a system of N spin-^ particles in a state thatis suitably analogous to dD, i.e. for which all 
of the terms have equal coefficients and for which the associated spin measurements 
admit an appropriate generalization of the Bernstein distribution. Having constructed 
the N-particle analogues of dD, we turn our attention to the broader class of N-particle 
Bernstein states for which the phase degrees of freedom are included. 


Fragile entanglement statistics 


6 


The generalization of Section | 2 ] to the four-particle case is not entirely obvious, 
since for N = 4 we wish to construct a four-particle entangled state and an appropriate 
associated set of observables such that the resulting probability distribution exhibits 
not only pairwise independence but also triplet-wise independence. Various 
intermediate cases can also be considered. For example, in the case of a 
four-particle system, one might consider the construction of a state for which 
pairwise measurements were independent, but triplet-wise and quadruplet-wise 
measurements were not. In what follows, we consider the situation where correlations 
emerge only at the highest level of participation. 

The fragility property likewise can be understood as holding in a suitably 
hierarchical sense. For N = 4 , the highest degree of fragility arises in the case of a 
four-particle state with the property that the associated three-particle reduced states 
are separable. But one can also consider the construction of entangled states with the 
property that the associated three-particle reduced states are not separable, and such 
that entanglement is destroyed only at the level of the two-particle reduced states. 
The entanglement of such states is thus in some sense more robust. In what follows, 
we shall be mainly concerned with fragility of the highest degree. 

An example of a four-particle Bernstein state that can be regarded as suitable 

generalization of the three-particle state is given as follows: 

= ^(iTTU) + imt) + inn) + inn) 

+ inn) + inn) + inn) + inn)). (s) 

It can be verified that if a Oz measurement is made on any one of the four particles, 
then the probabilities of the 'up' outcomes are 

1 

P(+ • • •) = p(« -!-••) = P(« •-!-•) = P(« ••-!-) = -. (9) 

Then we find that for the probabilities of two of the spins being 'up' we have 

1 

P(-|- -!-••) = P(-|- • + •) = P(+ • • +) = P(« + + •) = P(« + , +) = p(, , + ( 10 ) 

with similar results for the various other possible outcome pairs, thus allowing us 
to conclude that the outcomes of the spin measurements are pairwise independent. 
Likewise, for the probabilities of three of the spins being 'up' we find that 

1 

P(-i- -I- -I- •) = P(-i- -I- • -I-) = P(-i- • + +) = p(« + + +) = (11) 

8 

with similar results for the various other possible outcome triples, thus showing that 
the outcomes of the spin measurements are triplet-wise independent. Yet, when 
we make Oz measurements on all four particles we find that the outcomes are not 
independent. In particular, we have 

1 

P(-i- -I- -I- -I-) = 0 ^ P(-i- • • •) X P(« -I- • •) X P(« • -I- •) X P(« ••+) = —. 

16 


( 12 ) 
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One sees, therefore, that when a quantum system is prepared in the state 
the outcomes of measurements of the four spin operators cti^, a2z, Usz/ and 04^ 
have Bernstein statistics in an appropriately generalized sense. One can check, 
moreover, that the various three-particle reduced density matrices, for example, 
P234 = are separable—the proof for general N is given below—and 

hence that the state |Bd)) defined by dU) has the fragility property 

The combinatorial considerations leading to the explicit form of the fragile 
Bernstein state and, more generally, to the construction of a special Bernstein 
state |B^^^) for any number of particles, can be summarized as follows. 

First, since the probability of an 'up' outcome from a spin measurement has to 
be for each particle the |t) and |J,) states must appear with an equal number. Here 
we assume, for now, that the squared amplitudes of the coefficients appearing in the 
expansion of the state are all equal, as in ^ and dS]). Later we shall prove that this is 
a necessary condition for a Bernstein state. The number of terms in the expansion of 
|B^^^) must therefore be even, say, 2k for some k. 

Second, since each pairwise 'up' outcome has to have the probability |, the 
number of terms containing the pairwise-up combination I-- - t ••• t •••) has to be 
jk. Third, since the triplet positive outcome occurs with probability g, the number of 
terms of the form || | ^ has to be 

Continuing with this line of logic, we observe that the number of terms containing 
n spin-up states must be (and hence for n = N - 1 we deduce that the number 
of terms containing N - 1 spin-up states is {^)^~^k. But since the numbers of terms 
jk, jk, ■ ■ •, {^)^~^k, have to be integers, we deduce that k can only be or 2^“^. The 
latter case corresponds to an equal superposition of all the 2^ basis states for which 
the probability of N spins pointing in the up direction is 2 ~^, not zero, which has to 
be ruled out; so we deduce that k = 2 ^~^. This argument leads to the construction of 
the fragile Bernstein state |B^^^) for each N > 3 . 

As an example, let us show how one arrives at the state dll)- We know that 
is made up of 2fc = 8 terms. Four of them have the form |t • • •) and four of them 
have the the form !],•••). Of the first four, there are two of the form |tT • •) and 
the remaining two are of the form • •). One state in the former pair looks like 
ItTT •) and therefore must be ITTTi) since the particles carmot all be found in up 
states. The other state in that pair has the form Itti •) and is therefore ITTi-T) since the 
configuration with the first, second and fourth particles in up states is present exactly 
once. Now consider the two states of the form |tJ. • •). We see that one of them is 
of the type ItiT •) and the other one is of the type ItiJ. •), which follows since there 
must be exactly two terms with the first and third particles in up states. The first of 
these types can be completed as Iti-TT)/ since the configuration with first, third and 
fourth particles in up states has to be present once. But this implies that the other 
type, Iti-i. •), can be completed as Itslii) since there can only be two terms that have 
the first and fourth particle in up states. The correct expression for the remaining four 
terms, of type |J, • • •), can be identified in a similar maimer. This leads to dll) if one 
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takes the phase factors multiplying the various terms just identified each to be unity. 
The same kind of reasoning can be employed for general N. In the case of N = 5, for 
instance, one finds that the special Bernstein state takes the following form: 


|B(5)) = -(intti) + ITTTiT) + ITTUT) + ITiTTT) 

+ UTTTT) + imU) + intU) + ITUU) 

+ ITsLslslT) + UTTsli-) + IsLTslTsL) + UTsLiT) 

+1 i-i-TTi-) +1 slsLTslT) +1 i-slslTT) +1 ]•]•]•].].)) • 


(13) 


The fact that exhibits the required Bernstein statistics when Uz-up measurements 
are performed on each particle is ensured by the construction. 


4. Relation to GHZ states 

The special Bernstein state thus constructed for each N has the property that 
whenever one of the particles is traced out the resulting density matrix is separable. 
How to verify the separability of the reduced density matrices is not immediately 
obvious on account of the lack of a general sufficient condition for separability. It 
turns out, nevertheless, that these states can be turned into an N-particle GHZ state 
by means of a suitable change of basis. In particular, let us express \B^^) in terms of 
the basis by writing for each spin IT) = -^ (I ^) + I ^)) and |T) = -^ (I ^) - I ^))- 
Then we can show that can be expressed in the form 

|B™> = (14) 

We can show then, after a short calculation, that if one traces out any one of the 
particles, the resulting (N - l)-particle reduced density matrix is separable, which 
establishes the claim for general N. 

To show that the N-particle GHZ^ state ((14)1 emerges from the special Bernstein 
states, let us substitute | <—) = (|T) -i- |i)) and | -^) = (|T) - IT)) in (1X4)) . It is evident, 

due to the minus sign in ((14)) . that \B^^l) can be expressed as a linear superposition 
of all the 2^“^ states containing odd numbers of Uz-down states. With this in mind, 
we can ask what the probability is for a given particle being in the Uz-up state. One 
only needs to count the number of terms in the expansion thus obtained for \B^^) 
that have the given particle in the o^-up state. But the number of terms for which 
an odd number of the remaining N - 1 particles point in the cr^-down direction is 
2 N -2 follows that the probability of observing any given particle in the Uz-up state 
is as required. Similarly, the probability of observing any given pair of particles 
in Oz-vip states can be obtained by noting that each relevant term in the expansion 
of \B^^) corresponds to a partition of the remaining N - 2 particles into subsets of 
odd elements. There are such partitions, so the corresponding probability is |. 
Continuing in the same vein, one sees that the probability that any k < N - 1 given 
particles are in the Oz-np state is 2“^. However, not all of the N particles can occupy 
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az-up states simultaneously, since the term proportional to |tT • • • T) is absent from the 
superposition making up This establishes that, starting from the GHZx state 

dm), if we perform the reverse change of basis states for each spin, then the resulting 
state is the state that we obtained from the construction preceding (IT^ . 

The foregoing analysis leads to a comment on the following rather 
counterintuitive property of GHZ states. Suppose that we prepare an N-particle 
GHZz state, given by 



(15) 


Then we enquire about the statistics arising from a projective measurement of 
the spin for each particle in the x direction. Since irrespective of whether the spin 
is in the z-up state or the z-down state, the probability of observing the x-up state is 
equal, one might expect that the probability of an affirmative outcome should be j. 
This intuition indeed is correct. Next we ask what statistics arise for joint projective 
measurements of this type for pairs of particles. By continuation of the intuitive 
argument, and taking into account the fact that the state | GHZ^) is totally symmetric, 
we conclude intuitively that the probability of an affirmative outcome should be 
and this indeed turns out to be correct. However, when this logic is extended all the 
way to the N particle statistics, suddenly we find that intuition fails: the probability 
of an affirmative outcome is null, rather than 2~^. This probably comes as a surprise. 
Of course, if one were to perform the calculation by changing the basis in accordance 
with the above prescription, then one obtains a state of the form from which 

the correct statistics can be deduced. However, if one does not perform this analysis, 
the conclusion is not perhaps an obvious one. 

5. Construction of general Bernstein states 

The special Bernstein states N = 3,4,..., considered in the previous sections 

have been obtained under the assumption that the expansion coefficients are equal. 
What happens if we relax this condition? We know already that the inclusion of a 
phase factor in front of each of the basis states does not affect the Bernstein distribution. 
This follows from the fact that the basis states are orthogonal to one another, and that 
when the probabilities are calculated the phases cancel term by term. What is perhaps 
less obvious is that the squared amplitudes have to be the same. We proceed to prove 
that this is the case, and then turn to a discussion of the phase factors. 

To show that the magnitudes must be the same, we can use an inductive 
argument. Recall that for a system of N spin-^ particles there are 2^ basis states 
that can be constructed from the single-particle eigenstates. For the Bernstein 
distribution to hold, the all-up basis state | tt • • • T) rnust have coefficient zero. What 
about the basis states with one spin pointing down—for example, | tit • • • T) ? Glearly, 
there are N such states, one for each particle. The Bernstein distribution requires that 
the probability of finding any combination of N -1 spins in the up direction should be 
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2 i-N Hence, the coefficients of each of these N terms (up to phases) must be equal and 
given by Next, the contribution from these N terms towards the probability 

of finding N -2 spins pointing in the up direction is clearly 2 X 2^“^ = since for 

each fixed N - 2 spins pointing in the up direction there are two possibilities for the 
remaining two spins (for example, | tit ‘ ’ T) and | i,tT ‘‘ T) when we fix the last N-2 
spins). But for the Bernstein distribution to hold we require that the probability of 
finding N -2 spins pointing in the up direction must be 2^“^, which we already have 
achieved from the N terms, so there carmot be any other term containing N -2 spins 
in the up direction. It follows that the coefficients in front of the states containing two 
spin-down states (for example. Hit • • • T)) ruust be zero. Continuing with the analysis 
of these N terms, their contribution to the probability of finding N -3 particles in the 
up direction is clearly 3 X 2 ^“^ (for example ITTiT‘"T)/ITiTT’‘’T)/ and | ittt • • • T) 
all contribute towards the last N - 3 spins pointing in the up direction). The Bernstein 
distribution requires this probability to be 2^“^, so now we are short of terms. This 
can only be fixed by adding, for each fixed N - 3 spins in the up direction, a term for 
which the remaining three spins are pointing in the down direction (in the previous 
example, we need to add | iii'l ■ ■ ■ t))- The coefficients of these three-spin-down 
states must all, up to phases, be given by since 2^~^ - 3 x 2^“^ = 2^“^. 

We have shown that the coefficient of each term with a single spin pointing down 
or three spins pointing down is whereas all terms with no spins pointing 

down or two spins pointing down must be absent. The contribution from these 
existing terms towards the probability of finding N - 4 spins pointing up is given by 
4 X 2^“^ -I- 4 X 2^“^ = 2^“^, which is the number required by the Bernstein distribution. 
It follows that the terms with four spins pointing down must be absent. Continuing 
inductively, one sees that all terms containing even numbers of spins pointing down 
must be absent, and all terms containing odd numbers of spins pointing down must 
be present, with equal coefficients, up to phases. Hence, the construction of | ) 

presented in Section |3] gives the general form of a Bernstein state, up to phases. 

6. Manifold of general Bernstein states 

The foregoing analysis shows that the manifold Sn of projective N-particle Bernstein 
states is a torus in the projective Hilbert space PH^ of pure N-particle states. 

More precisely, since there are 2^“^ terms in the basis expansion of a Bernstein state 
in the Hilbert space of N spin-^ particles, there are 2^“^ phase factors that do not 
affect the Bernstein probability distribution. The associated 2^“^ - 1 relative phases 
parametrize a torus in PH . An interesting question to ask is whether local unitary 
transformations of the pure phase type ('local phase transformations') are transitive 
on !Bn- Equivalently, starting with a general Bernstein state, can the phase factors 
be eliminated by the application of a local phase transformation in such a way as to 
transform it into a special Bernstein state of the type constructed in Section with 
equal real coefficients? The answer turns out to be no, except in the case N = 3. 
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This can be seen as follows. We fix a point p in and apply a local phase 
transformation to it to obtain a new point Tp in Sn, which we call the transformed 
point. Thus, in each term of the expansion of the state vector representing the point 
p, we apply the phase factor exp{iak) if the state of particle k is |t), and we apply the 
phase factor exp{i(ik) if the state of particle k is \i). The components of the transformed 
state vector can be regarded as a set of homogeneous coordinates for Tp. In each term 
of the basis expansion of the transformed state vector, the phase of the corresponding 
term in the original state vector is shifted by a linear combination of the various ak 
and (ik- We call this combination the 'phase shift' of that term of the basis expansion. 

In the calculations below, it is convenient to remove a certain overall phase factor 
from the state vector in a way that ensures that the remaining coefficients depend on 
the various ak and only through the N phase differences defined by 6k = au- ^k for 
k = 1,2,.. .,N. For odd N we remove an overall phase factor of the form exp (Z,k=i (^k)- 
For even N, we remove an overall phase factor of the form exp (a^ + ll,k=i Pk)- In 
both cases, once the relevant overall phase factor has been removed, we redefine 
the homogeneous coordinates of the transformed point Tp accordingly The new 
homogeneous coordinates for Tp clearly belong to the same projective equivalence 
class as the original homogeneous coordinates for Tp. It follows that the submanifold 
Sjv(p) c Sn of projective Bernstein states that are attainable, starting at a given point 
p G Sn, by the action of local phase transformations is parametrized by a family of N 
phase differences, namely the 6k defined above. 

We shall show that for each choice of p there is a one-to-one invertible map from 
the points of an N-torus to the points of Sj^{p). Consider the N terms in the expansion 
of the transformed Bernstein state that contain a single spin-down particle. Denoting 
the phase shifts appearing in these terms by (pk where k = 1,2,...,N, we find that the 
various (pk can be expressed in terms of the various 6k by the relation 
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when N is odd, since for each term with one of the spins pointing down, the phase 
shift is given by a sum of all but one of the 6k. Similarly, when N is even we obtain 
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Note that the matrix of coefficients in (IT7|) is obtained by subtracting unity from each 
entry of the last column of the matrix of coefficients in (fT^ . For both N even and N 
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odd, the relevant matrix of coefficients is invertible, and one finds that 


61 ^ 

/ 

62 

1 

5 n-i 

1 

1 

1 — ' 


if N is odd, and that 
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(19) 


if N is even. Thus, the various 5k can be expressed in terms of the various (pk- 

The point p is left invariant under the given local phase transformation if and 
only if the phase shifts are integral multiples of 2n. More generally, two distinct local 
unitary transformations of the pure phase type transform p to the same point if and 
only if the associated (pk differ by integral multiples of 2n. This allows us to identify 
a toric periodicity structure in the space of 5k, which can be read off from the column 
vectors of the matrix in (fTS)) when N is odd and from the column vectors of the matrix 
dl^ when N is even. For instance, in the case N = 3 we see that (ITSb implies that 
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It follows that p is invariant under precisely those transformations that satisfy 
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for some collection of integers I, m, n. More generally, we see that 5k and 5k transform 
p to the same point if and only if 
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Once one has made the appropriate identifications, the points of the resulting torus 
are in one-to-one correspondence with the states that can be reached from the given 
Bernstein state p by use of local phase transformations. 
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In the case N = 4 we have 
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for some collection of integers k, I, m, and n. More generally, 6 k and 6 k transform p to 
the same point if and only if 
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Equation ([25l allows one to identify the periods of a torus, the points of which are in 
one-to-one correspondence with the states attainable from the given Bernstein state p 
by use of local phase transformations. 

More generally, the submanifold S^ip) of Bernstein states that can be reached 
from an initial state p by local phase transformations takes the form of a nonsingular 
embedding of an N-torus in with a periodicity structure similar to (|22|) . (|25| . For 
N > 3, this N-torus is of lower dimension than the Bernstein manifold Sn, which is a 
torus of dimension 2^“^ - 1. Thus, for N > 3 the phase factors cannot be eliminated 
by local phase transformations. But for N = 3 we have 2^“^ -1 = 3. 

An implication of this observation is that the states, together with those 

states accessible from a CFFL^ state via local phase transformations, form a proper 
subset of the space of general Bernstein states—the only exception being the case 
N = 3, for which GFFL^ states and general Bernstein states are equivalent modulo 
local phase transformations. This subset is of particular interest since all its elements 
satisfy the fragility property, namely, if one traces out the degree of freedom associated 
with any one of the particles, the resulting density matrix is separable. 


7. Inhomogeneous Bernstein states 

In the formulation of the quantum Bernstein distribution that we have considered so 
far, we have required that the relevant probabilities should be symmetric in the sense 
that the probability of any one given particle being in a spin-up state should be the 
probability of any two given particles being in up states should be and so on, all 
the way up to N - 1 particles. The resulting states are unique up to phase factors. 
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If we relax the symmetry property, and merely require that the measurement 
outcomes for any choice of 2,3,... ,N - 1 particles are independent, and that the 
measurement outcomes for all N particles are not independent, then there are many 
more possibilities. We shall refer to such states as an inhomogeneous Bernstein states. 
As an example, consider the state 

|Bf> = 9(lTU> + inT> + UTT>) 

+ V'!(l - 2<?) (l TU> +1 iTi) +1 UT>) + Vl-Sia-l) I ( 26 ) 

for 0 < q < For q = \ we have \B f For general q we find that 

2 

P(+ • •) = P(« + ^) = T{^^+) = q. (27) 

for measurements on a single particle, and that 

P(+ + •) = P(+ • +) = P(« + +) = q^ (28) 

for measurements on pairs of particles. Flence the outcomes are statistically 
independent. Yet, for measurements on three particles we have 

P(+ + +) = 0 ^ P(+ • •) X P(« + •) X P(« • +) = q^. (29) 

Thus we have a generalization of the Bernstein statistics. 

On the other hand, the reduced density matrix that results if we 

trace out the first particle is not separable unless q = \, and the same holds for the 
other two reduced density matrices. Flence the inhomogeneous Bernstein states do 
not inherit the fragility property. 


8. Remarks on Mermin's paradox 


The analysis that we have undertaken in the previous sections can be used to gain 
some insight into aspects of the so-called Mermin paradox [H. In particular, with 
a slight modification of our earlier calculations one can obtain a generalization of 
Mermin's result to any number of particles. 

It will be convenient to begin by reviewing the ideas of [H in a version presented 
by Peres | fT4| . A standard three-particle GFIZz state of the form (fT^ is prepared, and 
one considers the outcomes of product spin measurements made using four different 
combinations of directions, given by (x,y,y), (y,x,y), (y,y,x), and (x,x,x). The observable 
corresponding to the (x,y,y) combination, for instance, is given by Oix 0 02 y 0 o^y. The 
GHZ state is a simultaneous eigenstate of these four product observables. To see that 
this is the case, we recall that 
1 . 




1 





Fragile entanglement statistics 


15 


where we write \/') and |^) for the y-up and the y-down states. To calculate the 
eigenvalue of the GHZz state for the (x,y,y) observable, one expands the state of the 
first particle in the x-basis and the states of the remaining particles in the y-basis. 
Expressed in the new basis states, only terms with an even number of downward 
spins remain. Hence, the eigenvalue of <2) cr 2 y <2) o^y is 1. Analogous calculations 
for (y,x,y) and (y,y,x) show that the relevant eigenvalues take the value 1. To compute 
the eigenvalue for (x,x,x), one expands all the particles in the x basis. Only terms 
with an odd number of downward spins appear, leading to an eigenvalue of -1. A 
paradox then arises if one assumes that preassigned outcomes can be assigned to the 
various individual spin measurements, depending on which direction is chosen for 
the measurement axis. If we write Mnx for the outcomes of the Onx measurements, and 
for the outcomes of the Ony measurements, then it follows that mixm 2 ym^y = 1, 
miym 2 xm 3 y = 1, miym 2 ymi,x = 1, and mixm 2 xm 2 ,x = -1. On account of the fact that 
= m^y = 1, the product of the left sides of these four equations is +1, whereas the 
product of the right sides is -1, and hence we have a contradiction. 

In the calculations of the eigenvalues above we made use of the method of 
particle-wise coordinate rotation, examples of which we have already encountered. 
For instance, dl]), dH]) and (IT3l) are basis-transformed expressions of the GFFZx states 
for three, four and five particles, respectively. It turns out that an adaptation of the 
N-particle version of this method leads to a generalization of jU for N > 3. Consider 
an observable consisting of 2k factors of Oy and N - 2k factors of Ox- By expressing 
the respective 2k spin vectors in y-basis and the remaining N -2k spin vectors in the 
x-basis, one observes two patterns. If k is even, then only terms with an overall odd 
number of down spins remain. If k is odd, then all terms contain an even number of 
down spins. As before, this is inconsistent with preassigned values. 

For example, if N = 4 then two of the relevant observables are given by the 
combinations (x,x,x,x) and (y,y,y,y)- Here, k = 0 and k = 2 respectively, hence an odd 
overall number of downward spins remain, corresponding to an eigenvalue of -1. 
There are six additional observables, each with two factors of Oy. For these we have 
k = 1, and hence the eigenvalue is 1. We therefore obtain the relations 

mixmxmxmx = -1/ miymiym^ym^y = -1 (32) 

and 

mix^x^ymy = 1 , mixmym^xmy = 1 , mix^ym^ym^x = 1 , 

miym^xm^x^iy = 1, myrn^xm^yrn^x = 1, miym2ym^xmx = 1- (33) 

Unlike in the N = 3 case, where one encounters a single contradiction, now we 
encounter several distinct quadruplets of equations leading to contradictions. One 
such quadruplet, for instance, consists of the first equation of (l32l) , with equations one, 
two and three of (l33b . One can replace equations two and three by equations four and 
five to obtain a distinct quadruplet. In total, one finds eight distinct quadruplets. 

As another example, let us consider the case N = 5. The GHZ^ state is a 
simultaneous eigenstate of 16 observables containing even factors of Oy, namely, 1 
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of the form (x,x,x,x,x), 5 of the form (x,y,y,y,y), and 10 of the form (x,x,x,y,y), the last 
set of 10 having eigenvalue 1 and all others -1. That is, 

mixmxmxmixmsx = -1 (34) 


for the zero-y measurements. 


mixmiym^ym^ymsy = - 1 , miym2xm^ym^ymsy = - 1 , miym2ym^:cmymsy = - 1 , 

niiym2ym3ym^m^y = -1, miym2ym2,ym^ym^^ = -1 (35) 

for the four-y measurements, and 

mixmixm^xmymsy = 1 , mi^m2xm2,ym^^msy = 1 , mi^m2ym3^m4^m5y = 1 , 

miyniixmxmxnisy = 1 , mixniixmymymsx = 1 , mixniiymxmynisx = 1 , 

niiyniixmxmynisx = 1 , rnixmymymjixmsx = 1 , niiyniixmymAxmsx = 1 , 

niiyni2ymxm^rm5:c = 1 (36) 


for the two-y measurements. Now we encounter many contradictions. For each 
equation in (l34l) or (l35l) one can find one or more combinations of three equations in 
(|36ll that lead to inconsistencies, and likewise for each equation in (1^ one can find 
a combination of three equations in (IMll or (l35ll that leads to an inconsistency The 
number of such contradictions grows as N is increased. 

It remains a challenge to find practical applications of the unusual properties 
of Bernstein states. Since some of the configurations that we have considered give 
rise to situations where a fully entangled N-particle state separates if any one of the 
particles is traced over, the results may find applications in quantum cryptography, 
for example, in the area of 'secret sharing' protocols IITSl . An interesting application 
of the Bernstein distribution in coimection with the so-called 'principle of common 
cause' has been discussed in IIT^ . The approach to experimental realizations of 
Bemstein-type statistics in the context of photon polarization described in 
offers a further possible direction for practical developments. 
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